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! Abstract. We review the concept of finite-temperature form factor that was introduced re- 



cently by the author in the context of the Majorana theory. Finite-temperature form factors 
can be used to obtain spectral decompositions of finite-temperature correlation functions 
in a way that mimics the form-factor expansion of the zero temperature case. We develop 
the concept in the general factorised scattering set-up of integrable quantum field theory, 
list certain expected properties and present the full construction in the case of the massive 
Majorana theory, including how it can be applied to the calculation of correlation functions 
in the quantum Ising model. In particular, we include the "twisted construction" , which 
was not developed before and which is essential for the application to the quantum Ising 
model. 

> 
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^ ■ 1 Introduction 



Relativistic quantum field theory (QFT) at finite temperature is a subject of great interest 
which has been studied from many viewpoints (see, for instance, [1]). An important task when 
. studying a model of QFT is the calculation of correlation functions of local fields, which are 

related to local observables of the underlying physical model. For instance, two-point correlation 
functions are related to response functions, which can be measured and which provide precise 
information about the dynamics of the physical system at thermodynamic equilibrium. Although 
5^ i applications to particle physics often can be taken to be at zero temperature, many applications 

to condensed matter require the knowledge of the effect of a non-zero temperature on correlation 
functions. 

In this article, we review and develop further the ideas of [2] for studying finite-temperature 
correlation functions in integrable quantum field theory. 

In recent years, thanks to advances in experimental techniques allowing the identification 
and study of quasi-one-dimensional systems (see for instance O U] ) , there has been an increased 
interest in calculating correlation functions in l-|-l-dimensional integrable models of QFT (for 
applications of integrable models to condensed matter systems, see for instance the recent re- 
view [5]). Integrable models are of particular interest, because in many cases, the spectrum of 
the Hamiltonian in the quantization on the line is known exactly (that is, the set of particle 
types and their masses) , and most importantly, matrix elements of local fields in eigenstates of 
the Hamiltonian, or form factors, can be evaluated exactly by solving an appropriate Riemann- 
Hilbert problem in the rapidity space [H [71 [HI [9l [10] . 



*This paper is a contribution to the Proceedings of the O'Raifeartaigh Symposium on Non-Perturbative and 
Symmetry Methods in Field Theory (June 22-24, 2006, Budapest, Hungary). The full collection is available at 
|http://www.emis.de/journals/SIGMA/LQR2006.html, 
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At zero temperature, correlation functions are vacuum expectation values in the Hilbert space 
of quantization on the line. The knowledge of the spectrum on the line and the matrix elements 
of local fields then provides a useful expansion of two-point functions at space-like distances, 
using the resolution of the identity in terms of a basis of common eigenstates of the momentum 
operator and of the Hamiltonian. This is a useful representation because it is a large-distance 
expansion, which is hardly accessible by perturbation theory, and which is often the region 
of interest in condensed matter applications. Form factor expansions in integrable models at 
zero temperature have proven to provide a good numerical accuracy for evaluating correlation 
functions in a wide range of energies, and combined with conformal perturbation theory give 
correlation functions at all energy scales (an early work on this is 

One would like to have such an efficient method for correlation functions at finite (non-zero) 
temperature as well. Two natural (mathematically sound) ways present themselves: 

• "Form factor" expansion in the quantization on the circle. It is a general result of QFT 
at finite temperature [12^ [T3t [H] that correlation functions, at space-like distances, can 
be evaluated by calculating correlation functions of the same model in space-time with 
Euclidean (flat) metric and with the geometry of a cylinder, the "imaginary time" wrap- 
ping around the cylinder whose circumference is the inverse temperature. In this picture, 
one can quantize on the circle (that is, taking space as being the circle, and Euclidean 
time the line), and correlation functions become vacuum expectation values in the Hilbert 
space of this quantization scheme. Then, one can insert a resolution of the identity in 
terms of a complete set of states that diagonalise both the generator of time transla- 
tions and of space translations, as before, and one obtains a large-distance expansion for 
finite-temperature correlation functions. 

Unfortunately, the two ingredients required (the energy levels in the quantization on the 
circle and the matrix elements of local fields) are not known in general in integrable 
quantum field theory. We should mention, though, that exact methods exist to obtain non- 
linear integral equations that define the energy levels (from thermodynamic Bethe ansatz 
techniques, from calculations a la Destri-de Vega and from the so-called BLZ program), 
and that matrix elements of local fields were studied, for instance, in [151 [T6| [TTl I18j . 
Also, in the Majorana theory, the spectrum is known (since this is a free theory), and 
matrix elements of the primary "interacting" twist fields were calculated in [19^ [20] from 
the lattice Ising model, and in a simpler way in [21] directly in the Majorana theory using 
the free-fermion equations of motion and the "doubling trick" . 

• Spectral decomposition on the space of "finite-temperature states" . The concept of finite- 
temperature states, interpreted as particle and hole excitations above a "thermal vacuum" , 
was initially proposed more than thirty years ago and developed into a mature theory under 
the name of thermo-field dynamics [22l [231 [2l] i^^r a review, see for instance [25]). Ideas 
of H-l-dimensional integrable quantum field theory were not applied to this theory until 
recently. In [26], the concept of bosonization in thermo-field dynamics was studied, and, of 
most interest to the present review, in [2] the concept of finite-temperature form factor was 
developed - matrix elements of local fields on the finite-temperature Hilbert space. There, 
it was studied in depth in the free Majorana theory, both for general free fields (finite 
normal-ordered products of the free Majorana fermion fields - including the energy field) 
and for twist fields. It was found that a Riemann-Hilbert problem again characterises 
finite-temperature form factors of twist fields, but that this Riemann-Hilbert problem 
presents important modifications with respect to the zero-temperature case. Solutions 
were written explicitly for primary "order" and "disorder" twist fields, and the full finite- 
temperature form factor expansions of two-point functions were written and interpreted 
as Fredholm determinants. 
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An interesting discovery of [2] is that these two methods are actuahy related: it is possible 
to evaluate form factors on the circle from (analytical continuations of) the finite-temperature 
form factors, and the analytical structure of finite-temperature form factors (and of the measure 
involved in the expansion of correlation functions) is directly related to the spectrum in the 
quantization on the circle. This provided a new way of evaluating form factors of twist fields 
on the circle, and most importantly, gave a clear prescription for the integration contours in the 
finite-temperature form factor expansion (naively plagued with singularities) . The requirements 
brought on finite-temperature form factors by this relation constitute, in a way, a generalisation 
of the modularity requirements found in conformal field theory for constructing correlation 
functions from conformal blocks. 

It is important to realise, though, that both expansions for correlation functions are not 
equivalent. The first one gives an expansion at large (space-like) distances, whereas the sec- 
ond can be used to obtain both large-distance and, expectedly with more work, large-time 
expansions. Indeed, the finite-temperature form factor expansion can naturally be deformed 
into an expansion in the quantization on the circle through the relation mentioned above [2]. 
It is expected that it can also be manipulated to obtain large-time behaviours. A manipu- 
lation of this type was done in [27j. There, going in reverse direction as what is described 
in [2], the expansion on the circle in the quantum Ising model was first deformed into a kind 
of finite-temperature form factor expansion (without being recognised as such), which was then 
used to obtain large-time dynamical correlation functions in a certain "semi-classical" regime 
(partly reproducing earlier results of [28] and [29]). This manipulation, however, neglected 
contributions that may change the asymptotic behaviour, and a more complete derivation of 
the large-time behaviours from finite-temperature form factor expansions is still missing. In 
particular, for the quantum Ising model, the Fredholm determinant representation of [2] and 
those obtained in the present paper may be of use, following the technology reviewed in [30] 
(work is in progress [31]). 

It is worth noting that the method we review here is not adapted to providing information 
about one-point functions at finite-temperature. Various works exist concerning such objects 
\32 \ 133 1 [M]. Work [32] is interesting in that it uses the knowledge of the zero-temperature form 
factors in order to deduce the finite-temperature one-point function of the energy field. The idea 
is to "perform" directly the finite-temperature trace from the known matrix elements. A regu- 
larisation is necessary, but the finite- volume one seems impossible to tackle. A certain convenient 
regularisation was proposed there and shown to reproduce the known finite-temperature average 
energy. The idea of using this regularisation for multi-point correlation functions has been 
suggested and we are aware of results in this direction [3^ , but it is not yet understood why in 
general this should work. 

Let us also mention that correlation functions of twist fields in the Major ana theory can 
be obtained as appropriate solutions to non-linear differential equations [36]. But at finite 
temperature, or on the geometry of the cylinder, these equations are partial differential equations 
in the coordinates on the cylinder [37^ [38] [39] , and do not immediately offer a very useful tool 
for numerically evaluating correlation functions, neither for analyzing their large-distance and 
large-time behaviours. 

The theory developed in [2] for the Majorana case is still incomplete. Twist fields present 
certain complexifications at finite temperature that are not present at zero temperature, and, in 
order to describe all correlation functions, one also needs a "twisting" of the construction of [2], 
as it was mentioned there. In addition, certain exponential pre-factors were omitted in [2]. 
These two aspects are in fact essential for applications of the results in the Majorana theory to 
correlation functions in the quantum Ising model. 

In this article we will review the ideas of [2], by developing them in the general factorised 
scattering context of integrable quantum field theory, and complete the work for the Majorana 
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theory. We will deduce many of the immediate properties that arise in the general context 
for finite-temperature form factors, drawing on the ideas of [2], and we will present both the 
untwisted and the twisted constructions. We will recall the results for the Majorana theory, 
and extend them to the twisted case, finally giving the explicit representation for correlation 
functions in the quantum Ising model. 

The article is organised as follows. In Section [2] we review the form factor program at zero 
temperature, and in Section [3] we recall basic results about finite-temperature correlation func- 
tions. Then, in Section^ we describe the concept of finite-temperature states using the language 
of factorised scattering in integrable QFT, we introduce the concept of finite-temperature form 
factor and we describe the resulting expansion of correlation function. We also present the 
ideas underlying the relation between finite-temperature form factors and matrix elements in 
the quantization on the circle, still in the general context. In Section \5\ we develop the basics 
of the twisted construction. In Section [6l we present certain formal results about the space of 
finite-temperature states, and in particular, we deduce a generalisation of the idea of "mapping 
to the cylinder" that one uses in conformal field theory in order to study correlation functions 
at finite temperature (again, this is a generalisation of ideas of [2]). Finally, in Section [Tj we 
recall and extend the results of [2] for the Majorana theory and its connection to the quantum 
Ising model. 

2 Review of the zero-temperature form factor program 
in integrable quantum field theory 

The Hilbert space of massive relativistic quantum field theory is completely specified by fixing 
the set E of particle types of the model. In 1-1-1 dimensions, every Hamiltonian eigenstate is 
then described by choosing k £N particle types and by associating to them k real numbers, the 
rapidities: 

1^1, • • • , dk)ai,...,ak 

with ai £ E and 9i S M (and the order of the rapidities/particle types is irrelevant - hence a basis 
is obtained by fixing an ordering of the rapidities). The Hamiltonian H and the momentum P 
act diagonally on these states. In order to fix their eigenvalues, one only has to fix the masses 
nia £ K"*" for every particle type a £ E. The eigenvalues are then 

k k 
H : Ek = y^^rug^ cosh 9 j, P : = ^ sinh 6*, . 

i=l i=l 

Other symmetries of the model also act diagonally, and their eigenvalues are fixed by choosing 
charges associated to the various particle types. 

There are many possible bases of the Hilbert space, all described as above. Two are of 
particular importance: the in basis and the out basis. They describe, respectively, particles of 
the given types and rapidities far in the past, and far in the future (in non-integrable models, one 
should really include the additional dependence on the impact parameters). The far past and 
the far future are regions in time where all particles are so far apart that they do not interact, 
and can be described as freely propagating. The overlap between the in basis and the out basis 
gives the scattering matrix: 

1^1, ^2, . . .)S2,... = E / ^^i'^^2 • • • Sai:ar:::i9i, 9^, 

II " 



; 612, . . .)|( 



''li^2) 



, (put) 
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where the number of particles in the in state and in the out states is generically different. The 
structure of the Hilbert space and the Hamiltonian describe the particles and their propagation, 
but it is the scattering matrix that encodes the interaction, and in particular, the locality of 
relativistic quantum field theory. 

In integrable quantum field theory, the scattering matrix can be determined from the physical 
requirements of unitarity and crossing symmetry, from the integrability requirement of factori- 
sation and the lack of particle production, and from minimality assumptions and the "nuclear 
democracy" (every pole has a physical explanation through resonances from particles already 
in the spectrum). All scattering processes can then be described using only the two-particle 
scattering matrix Salfa^i^i ~ ^2), Oi > 62'- 




bi,b2 



It is convenient for this purpose to introduce the Zamolodchikov-Faddeev algebra (from now on 
in this section, summation over repeated indices will be implied) 

Z"n^l)^"H^2) - ^feX^^l - ^2)^'H^2)^'n^l) = 0, 

ZaMZaM) - S'aliliOl " ^2)^6, (^2)^6, (^l) = 0, (2.1) 

z--{e,)z,,ie2) - s'^i';i{e2 - ei)z,,{e2)z'^H0i) = sz 5{e, - e^). 

The in basis and the out basis are then two bases for the same Fock space (actually, a general- 
isation of the concept of Fock space) over this algebra, defined simply by different ordering of 
the rapidities: 

Z"^ (0)1 vac) = 0, 

1^1, • • • , Ok)tl^a, = Za, {Ol) ■ ■ ■ ^a,(efc)|vac) (^i > • • • > ^fc), 

1^1, • • • , Ok)tla, = Za, (Ol) ■ ■ ■ Za, m\y^c) (^i < • • • < 0fc). 

The natural Hermitian structure on this space gives {Z"'{0))^ = Za{9). 

Once the Hilbert space has been identified with the Fock space over the Zamolodchikov- 
Faddeev algebra, the algebra elements Z°'{9) and Za{9) become operators with an action on the 
Hilbert space. It turns out, from expected properties of quantum field theory, that they induce 
very nice properties on the objects (form factors) 

<,...,a,(^i, ■■■,dk)^ (vac|O(0, 0)Za, i9i) ■ ■ ■ Z,J0fc)|vac), 

where 0{x, t) is a local field of the model. Indeed, these objects, defined here for real rapidities, 
actually are (by analytical continuation) meromorphic functions of the rapidities. They can be 
determined through a set of analyticity requirements and through the recursive determination 
of the residues at the poles (form factor equations) [8} [TO]: 

1. Meromorphicity: as functions of the variable 9i — 9j, for any i,j £ {I,... ,k}, they are 
analytic inside < lm(0j — 9j) < 27r except for simple poles; 

2. Relativistic invariance: 

Fa„...,aA(^l +P,...,9k+P)= e^(^)^i^f?,...,a,(^l, • • • 

where ^(0) is the spin of O; 
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3. Generalized Watson's theorem: 

^ai,...,aj,aj+i,...,aS^^^ ■ ' ' • • • ) ^fc) 

= Sa/XV'l i^j ~ ^j+i)^ai,...,bj+i,bj,...,ak (^1) • • • i %5 ■■■,0k)] 

4. Locality: 



where fo is 1 if O is fermionic, if it is bosonic, ^ is the fundamental field associated to 
the particle a^, and uj{0, ^) is the semi-locality index (or mutual locality index) of O with 
respect to ^ (it will be defined in Subsection 13. 3p : 

5. Kinematic pole: as function of the variable 9n, there are poles at Oj+in for j G {1, . . . , k—1}, 
with residue 



Faj_,...,aj,...,ak^i{Gl, ... ,6j,... , 9k-l) 



- (-i)/o/*e2--«^''f)5^ti • • • ^^s'i^;^~\(o, - e^^^^^' 

where a hat means omission of the argument, and Ca^,bj is the conjugation matrix. 

6. Bound-state poles: there are additional poles in the strip < lm(0j — Oj) < vr if bound 
states are present, and these are the only poles in that strip. 

Form factors can in turn be used to obtain a large-distance expansion of two-point correlation 
functions of local fields: 



(vac|Oi(x,t)O2(0,0)|vac) = f; J] M^i ^ " ^^fc ^-^^E, i^.+»^E,P. 

fc=Oai,---,afc 

X (vac|O(0, 0)1^1, ... , ek)a,,...,a, au...,a,{9i, ■■■, 0k\O{0, 0)|vac). 

A large-distance expansion is effectively obtained by shifting all rapidity variables by 7r/2 in 
the positive imaginary direction, and by using relativistic invariance. This gives a formula 
which looks as above, but with the replacement e~**^J '^■''^^^^j^^ ^ e 5Dj '"a^ co&h(ej ) ^j^g^.^ 
r = \/ — t'^. It turns out that this is numerically extremely efficient in most integrable models 
that were studied. 

3 Finite temperature correlation functions 
3.1 Traces 

Physical correlation functions at finite temperature are obtained by taking a statistical average 
of quantum averages, with Boltzmann weights e~^^ where E is the energy of the quantum state 
and P is the inverse temperature. They are then represented by traces over the Hilbert space: 
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Since all matrix elements of local fields are known in many integrable models, it would seem 
appropriate to write the trace as an explicit sum over all states of the Hilbert space, and to 
introduce resolutions of the identity between operators inside the trace, in order to evaluate 
finite-temperature correlation functions. However, this method does not account correctly for 
the fact that at finite temperature, states that contribute to the trace are very far from the 
vacuum. Yet, it turned out to give good results in the case of correlation functions with only 
one operator [331 EH [32] • 



3.2 Quantization on the circle 

On the other hand, traces as above can be represented by vacuum expectation values on the 
Hilbert space of quantization on the circle of circumference p. Indeed, a consequence of the 
imaginary-time formalism [12] is the Kubo-Martin-Schwinger (KMS) identity \13\ I14j. 

{{0{x,t) . . .))f, = - ■ ■ ■))(,, (3.2) 

where (— l)'^'' is a sign accounting for the statistics of O (it is negative for fermionic operators 
and positive for bosinic operators), and where the dots (• • • ) represent local fields (that are also 
local with respect to O) at time t and at positions different from x. Then, finite-temperature 
correlation functions can be written as 

((0(T,ix)---))^ = (e*-^/2...) ^(vac|a(x,r)---|vac);3, (3.3) 

where s is the spin of O, and there are factors for all operators in the correlation 

function. The operator ©(x, r) is the corresponding operator acting on the Hilbert space Hp of 
quantization on the circle, with space variable x (parameterizing the circle of circumference 0) 
and Euclidean time variable r (on the line). The vector |vac)^ is the vacuum in this Hilbert space. 
Below, we will mostly be interested in fermionic models, that is, models with a "fundamental" 
fermion field (which creates from the vacuum single-particle states). For such models, one 
can think of at least two sectors in the quantization on the circle: Neveu-Schwartz (NS) and 
Ramond (R), where the fundamental fermion fields are anti-periodic and periodic, respectively, 
around the circle. The trace (|3.ip with insertion of operators that are local with respect to the 
fermion fields naturally corresponds to the NS sector due to the KMS identity. This is the sector 
with the lowest vacuum energy. 

The representation (13. 3p immediately leads to a large-distance expansion of finite-temperature 
correlation functions, through insertion of the resolution of the identity on the Hilbert space W^: 

gE,".'^+(A£-B„i,...,„,)r 

;j(vaci|O(x,r)O(0,0)|vac2)/3 = ^ (3.4) 

A;=0 ni,...,nfc 

X /3(vaci|(!)(0, 0)|ni, . . . ,nfc)/3 /3(ni, . . . , nfc|O(0, 0)|vac2)/3, 

where the eigenstates of the momentum operator and of the Hamiltonian on the circle are 
parametrized by discrete variables ?ij's. The vacua |vaci)/3 and |vac2)/3 may be in different sec- 
tors, and these sectors may be different than the sector where the excited states |ni, . . . 
lie (this situation occurs when considering semi-local operators as is recalled in the Subsec- 
tion 13.31 below) . The quantity AiS is the difference between the vacuum energies of the vacuum 
state |vaci)/3 and of the vacuum above which the states |ni, . . . , nfc)/3 are constructed. The 
states |ni, . . . and the excitation energies -Eni,...,nfe may also depend on additional discrete 

parameters (quantum numbers, particle types), on which one has to sum as well. 

This form is valid for any integrable model on the circle. However, this Hilbert space has 
a very complicated structure, even in integrable quantum field theory; for instance the energy 
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levels -E'ni,...,nfc are not known in closed form. Also, there is yet no known procedure in general 
integrable quantum field theory for evaluating form factors on this Hilbert space. Moreover, 
this representation does not provide large-time (real time) expansions, since it inherently gives 
finite-temperature correlation functions in Euclidean time. 



3.3 Semi-locality: U{1) twist fields 

If the model we are considering has internal global symmetries, then there are local twist fields 
associated to them. Twist fields are of interest, because they usually correspond to some order 
parameter. We will clarify the correspondence between order /disorder parameters in the quan- 
tum Ising chain and twist fields in the Section 17.61 The first appearances of certain twist fields 
in the context of the Ising statistical model can be found in [401141]. but we are going to describe 
twist fields here in more general terms (see, for instance, the lecture notes |42j). 

Twist fields are not local with respect to the fundamental fields associated to a given particle 
(but are with respect to the energy density). If the symmetry to which they are associated is 
C/(l) or a subgroup of it (and if the fundamental field transform by multiplication by a phase), 
then the twist fields are said to be semi-local with respect to the fundamental field. In the 
quantization scheme on the line, a twist field, which we will generically denote by a, gives rise 
to a pair of operators, which we will denote by cr+(x,t) and cr_(x,t), having a cut towards 
the right (positive x direction) and towards the left (negative x direction), respectively. These 
operators lead to the same correlation functions at zero temperature. 

When considering correlation functions at finite temperature, things are more subtle. The 
exact shape of the cuts are unimportant, but it is important if the cut is towards the right or 
towards the left. This is because the insertion of an operator a±{x,t) that is semi-local with 
respect to the fundamental field ^{x,t) may affect the vacuum sector in the correspondence to 
expectation values in the quantization on the circle. Semi-locality can be expressed through the 
exchange relations 

^{x, t)a+{x', t) = (-l)^*^''e-2-*'^®(^"^'V+(x', t)^{x, t) {x / x') (3.5) 

and 

^{x,t)a.{x',t) = (-l)^*^-e2"^'^®(^'-^'V_(x',t)^(rE,t) {x / x'), (3.6) 

where @{x) is Heaviside's step function and io is the semi- locality index associated to the pair 
(^', cr). Taking here the fundamental field to be fermionic (because this is what will be of interest 
in the following - the case of bosonic fundamental fields is straightforward to work out), it is 
a simple matter to generalise the KMS identity to (using /ij, = 1) 



{{^ix,t)a+ix',t)---))fs 
{{^{x,t)a.{x',t)---))fs 



-e-^^'^{{^{x, t - i/3)a+{x', t) ■ ■ •))^ {x ^ oo), 
-((^'(x, t - i(3)a+{x', t) ■ ■ ■))f3 {x -oo) 

-((^'(x, t - i(5)a.{x', t) ■ ■ {x oo), 

{{^{x,t-if3)a-{x',t)---))p (x^ -oo). 



-e 



where the dots (• • • ) represent fields that are local with respect to the fermion field ^, at time t 
and at positions different from x. Then, in the correspondence of the trace with a vacuum 
expectation value in the quantization on the circle, one of the vacua will be in a different sector, 
in accordance with these quasi-periodicity relations. Denoting by Ivacj,)^ the vacuum in the 
quantization on the circle with quasi-periodicity condition ^' i-^ e^^'^*'^^' around the circle in the 
positive space (x) direction, we have 

((cT+(r,ix)---))^ = (e*"^/2...) /3(vaci+^|a(x,r)---|vacNs)/3 (3.7) 
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and 



((a_(r,ix)-..))^ = (e*-^/2...) ^(vacNs|^(x, r) 



|vaci_^)^, 



(3.8) 

represent operators that 



where |vacNs)/3 = I vac 1)^3 is the NS vacuum, and where the dots (• 
are local with respect to the fundamental fermion fields. 

With many insertions of semi-local operators, similar phenomena arise. This change of the 
vacuum sector has an important effect: under translation in the x direction, the insertion of 
an operator a±{x^t) inside a trace produces an additional real exponential factor, due to the 
difference between the vacuum energies of the different sectors; that is, the trace is not trans- 
lation invariant. It is convenient to represent this lack of translation invariance of traces, in the 
case where many semi-local operators are multiplied, by considering "modified" transformation 
properties of this product of semi-local operators. Consider the product 0"^^^ (xi, ti) • • • (T':^^{xk,tk) 
where r]i = i: and we have indicated explicitly the semi-locality indices u>i. Then, inside traces 
at temperature P with insertion of operators that are local with respect to the fundamental 
fermion field, we have 



AS = £ 



{xi,ti) • ■ ■ a'^^lixk.tk)- 



1 ^ 

i=l 


-8 


1 

2" E^^ 

i = l 


»li = + 




')i = - 



(3.9) 



where £[u\ is the energy of the vacuum |vac 



yip- 



4 A space of "finite-temperature states" 
in integrable quantum field theory 

In [2], it was suggested that the difficulties in obtaining large-distance or large-time expansions of 
finite-temperature correlation functions can be overcome by constructing a finite-temperature 
Hilbert space in terms of objects with nice analytic structure, in analogy with the zero-tem- 
perature case. The program was carried out explicitly in the free massive Majorana theory 
(considering, in particular, "interacting" twist fields). As we said in the introduction, the idea 
of a finite-temperature Hilbert space is far from new, but it is in [2] that it was first developed 
in the context of an integrable quantum field theory. 

4.1 General idea 

The idea of the construction is simply to consider the space C of endomorphisms of as a Hilbert 
space with inner product structure 



{A,B) 



Tr [e-P^A^B) 



Tr {e-P") ■ 

This Hilbert space is known as the Liouville space [l3]. Note that 
{A,Br = {B,A). 

There is then a (generically) one-to-two mapping from End('H) to End(£): to each operator C 
acting on TL, there are two operators, 4>l{C) and (j)R{C), acting on C, defined respectively by 
left action and by right action of C as follows: 

, , , Tr(e~/^"A^CB) , , , , Tr(e-^^A^BC) 
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In particular, if Q is a generator of a symmetry transformation on 7i, then — (f'RiQ) is 

the generator on C. The set of ah operators on C that are in the image of at least one (pi 
or (pji will be denoted EndLij(£). 

The main power of this construction, from our viewpoint, is the possibility to obtain large- 
distance or large-time expansions at finite temperature, in analogy with the zero-temperature 
case, using a resolution of the identity on the space C. Indeed, suppose we have a complete set 
of orthonormal operators D{9i, . . . , 9k), 9i > ■ ■ ■ > 9^ k G'N: 

{D{ei, 9k),Die[, 91)) = 5k,i 5{0i - e[) ■ ■ ■ 5(9^ - 9',). 

Then, we can decompose any inner product as a sum of products of inner products: 

oo „ 

{A,B) = Y, / d9i---d9k{A,D{9,,...,9k)){D{9i,...,9k),B). 

This is a non-trivial relation equating a trace on the left-hand side to a sum of products of traces 
on the right-hand side. 



4.2 A natural basis 

A natural complete set of operators in integrable quantum field theory can be obtained as follows 
(more precisely, one should consider an appropriate completion of the set below). First, define 
a larger set of particle types Ec = E ® E, the elements being the couples a = (a, e) for a & E 
and e = lb. For notational convenience, define = Za e = + and Z^ = Z"" if e = —. Then, 
we have a complete set with 

Dai,...,ak{Ol, ■ ■ ■ ,9k) = Za^{9i) ■ ■ ■ Za^{9k) 

for 



ai£Ec, 9i>--->9k£R, k£N. (4.1) 

In fact, it will be convenient to define the operators -Dqi,...,^^. (^i, • . . , 9k) for any ordering of the 
rapidities, and to define them as being exactly zero when two rapidities collide (in order to avoid 
overlap with operators with smaller k): 

{ Z^,{9,)---Z^^{9k), 9,^9, y 1^3, 

ID, 9i = 9j for some i f= j- 

These operators will form a very useful set if the matrix elements (with 0{x, t) a local operator 
and 0^{x,t) its Hermitian conjugate on Tl) 

{0\x,t),D^,_^^{9i,...,9k)) 

have simple analytical properties; for instance, if the value of this function for a certain ordering 
of the rapidities is the one obtained by analytical continuation from that for a different ordering. 
Then, it may be possible to write down equations similar to the zero-temperature form factor 
equations (what we will call the "finite-temperature form factor equations") and to solve them. 
Two clues suggest that this may be so, at least for models with diagonal scattering (see below): 
first, these objects specialise to the zero-temperature form factors when the temperature is sent 
to zero and all signs e^'s are set to -|-, and second, in [2] the finite-temperature form factor 
equations were indeed written and solved in the free massive Majorana theory (this will be 
recalled in Section [7]). 
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Although the operators (|4.2p form a complete set, they are not orthonormal. In the case of 
diagonal scattering (and this is the only case we will consider from now on) 

(without summation over repeated indices), it is possible to write down all inner products in a 
simple way, using the Zamolodchikov-Faddeev algebra (j2.ip and the cyclic property of the trace: 

iDai,...,a,iOi, ■ ■ ■ ,^fc),^a;,...,a;(^l, • • • , O'l)) = 5k,i n ^ _\_iJa-''^-uP^a, coi^i' (4-3) 

i=l ^ ' 

where from unitarity, (— l)-'^" = «Sa,a(0) = ±1 (this corresponds to the statistics of the particle of 
type a, as an asymptotically free particle). Here, we have assumed the ordering ()4.1I) for both 
members of the inner product. 

Note that simple "crossing" relations hold for operators in End/j/,(£): 

= e'^'^^ '°'''''(I?«„...,„„^;(ei, . . . , e,, 9[) , H{A)D^>^_^,je[, . . . , OU)) {e[ / Vz) 



and 



{Da„...,a, {Ol,..., 9k), MA)Da'„...,a[ (^1, • • • , (4.5) 



where (a,e) = (o, — e). 

4.3 Finite-temperature form factor expansion 

Inverting (j4.3p . and using the fact that the operators (j4.2p are eigenoperators of both the Hamil- 
tonian and the momentum operator, we get a spectral decomposition for two-point functions 
{{Oi{x,t)O2{0,0)))i3 (finite-temperature form factor expansion). In order to simplify the dis- 
cussion, we assume that x > 0. This can always be achieved by taking complex conjugation if 
necessary: {{Oi{x,t)O2{0,0)))*p = {{Ol{0,0)Ol{x,t))) (3 (a slightly different formulation holds 
with X < 0). We have 

{{Oi{x,t)O2{0,0)))f3 (4.6) 

7/1 in T-r l~/<2, , ei(ima xsinhdi —ima t cosh 9i) 

= ^^'^E E / 



kin 1 - (-l)^''^e" 
i=i ^ 

X F^l..,a, (^1, • • • , 0k; l3)F%__^^{ek, . . . , ^i; /?), 
where we have defined finite-temperature form factors as the normalised matrix element^: 

FZ,...,a,{Oi,...,ek;P) (4.7) 

n 

= n V^^^' (l - (-1)^"' ^-^.0m..^oshB,y^ (Ot(0,0) , (01, . . . ,0fc)). 

i=l 



^Note that by definition of the basis of states in the function -F^^ {9i, . . . , 9n\ P) has no delta-function 
contributions at coUiding rapidities. 
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This normalisation is for later convenience (one may call it the "free field" normalisation). It 
leads in particular to the following identity, which we have used to write the expansion: 



which essentially follows from (j6.2p below. In the expansion (j4.6p . we have symmetrised over 
the orderings of rapidities. The quantity A£ is non-zero whenever Oi is a twist field aj^j^, and 
is given by 



AS 



+ UJl 

-£ 



UJ-\ 

2 



im = +), 
ivi = -), 



where <5[z^] is the energy of the vacuum jvacj^)^ (see the discussion around (|3.9p ). 
Some comments are due: 

• As we said, the factor ^ is present whenever the operator Oi is semi- local, Oi = (y^, 
with respect to the fundamental fermion field. It is as in (13. 9p with k = 1, and with uji = uj 
and rji = Tj the semi- locality index and cut direction, respectively, of Oi. It occurs because 
the operators Za{6) can be expressed through integrals of the fundamental fermion field. 
To be more precise, in order to deduce it from the discussion around (j3.9p . one has to 
assume that although these integrals extend to ±00 in the x direction, they only produce 
excited states, without changing the sector. The presence of this exponential factor can in 

fact be shown for the finite-temperature form factors of the order field a^. in the Majorana 

1 

theory. Indeed, as we said, in the Majorana theory the traces -Fq{';...,q„ (^1, . . . ,9n', P) were 
shown in [2] to satisfy a set of recursive relations which ultimately relate them to the 
one-point function (in the case of the order field). Slightly generalising the derivation 
to include an x and t dependence, this accounts for the phase factors above. On the 
other hand, the one-point function of a twist field is not translation invariant, as is clear 
from (j3.7p and ()3.8p . the transformation property being as in (j3.9p . This is what accounts 
for the real exponential factor (this factor was missing in [2], because the one-point function 
was considered translation invariant). 

• When both Oi and O2 are semi-local with respect to the fundamental fermion fields, the 
finite-temperature form factor expansion (j4.6p is valid only when the cut of Oi extends 
towards the right (positive x direction) and that of O2 extends towards the left (negative 
X direction). This will be justified in the Subsection 14. 4[ Note that with this prescription 
on the directions of the cuts, one produces the correlation functions 



/3(vaci+^JOi(x,r)a2(0,0)|vaci 



when written in the quantization on the circle, with x = t and t = ix. This is a restriction, 
as not all vacua on Ti/^ can be obtained. For instance, one would like to evaluate correlation 
functions of twist fields with the NS vacuum. This restriction will be lifted in Section [5j 



Had we not put a small imaginary part to the rapidities in the integrals, the expansion (j4.6p 
would have been plagued by singularities: as the rapidity associated to a particle of type 
(a, e) becomes equal to that associated to a particle of type (a, — e), poles are expected to 
appear in the finite-temperature form factors (kinematic poles). This expectation comes 
from the intuition from zero-temperature form factors, and from the fact that these singu- 
larities indeed occur in the finite-temperature form factors of twist fields in the Majorana 
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theory, as was calculated in [2]. There it was shown that a proper solution is obtained by 
slightly deforming the integration contours as indicated above: 



That this is the right prescription still in the interacting case will be argued in Subsec- 
tion 1131 

• It is important to realise that the expansion (|4.6p is not directly a large-distance or a large- 
time expansion. But it can be made so as follows (for large-time expansions, this requires 
some work). First, with further displacement of the integration contours in the directions 
of (j4.8p . more precisely with lm(6j) = ej 7r/2, the expansion (j4.6p becomes an expansion at 
large — t^ (recall that we consider x > for simplicity). In order to perform this contour 
displacement, one needs to know about the analytical structure of the integrands; this will 
be briefly discussed in Subsection 14.41 Second, the integrals involved in ()4.6p can be made 
convergent in time-like regions — > by deforming the contours in the following way: 

\t\ > \x\, t > : Im(6'j) = -sgn{Re{0j)) ejO+ , 
|t| > t < : Im(%) = sgn(Re(6lj)) 6^0+. 

These deformations necessitate the addition of residues coming from the kinematic poles. 
These residues will lead to powers of the time variable, which will need to be re-summed. 
Note that it was assumed in |27j that considering the contributions near the singularities 
at colliding rapidities, the expansion gives the leading in some semi-classical region, which 
should include a large-time limit — — > oo. The full contour deformation should give 
a definite answer as to the large-time dynamics (work is in progress [31]). 

• When calculating the spectral density (from the Fourier transform of the two-point func- 
tion), the expansion (14. 6p does produce an expansion with terms of lesser and lesser impor- 
tance as the particle number is increased, at least for large enough energies. However, one 
does not have the situation where the spectral density is known exactly up to a certain 
energy depending on the number of particles considered, as happens at zero tempera- 
ture. It would be very interesting to have a full analysis of the spectral density at finite 
temperature. 

4.4 Prom finite-temperature states to the quantization on the circle 

A great part of the structure of the finite-temperature form factor expansion can be understood 
according to the following idea. 

Suppose that we have a model of quantum field theory; more precisely, let us consider a statis- 
tical field theory, on a space with Euclidean signature. Let us quantize it with a certain choice 
of space X of infinite extent, and Euclidean time . If we were starting from a Lorentzian 
quantum field theory, with real time t, we would just be considering the Wick rotated variable 
t = -it^. Then, the Hilbert space is the space of field configurations on x, with appropriate 
asymptotic conditions. On this space, we choose a vacuum |vac) such that correlation functions 
are vacuum expectation values. Now suppose that a basis of states is chosen such that the 
generator of x translations is diagonalised. The operator producing x translations is unitary, of 
the form e~^^^, where P is the Hermitian generator. Suppose that the states are parametrised 
by the real eigenvalues p of the operator P. Since space is of infinite extent, then p takes all 
real values. 

Then, formally, if we were to "analytically continue the theory" towards positive imaginary 
eigenvalues p = iE, the operator producing x translations would have the form e^^ for some 
Hermitian H with eigenvalues E. The claim is that the operator H is still the generator of x 



Im(%) = ejO 
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translation, but now in a different quantization scheme (that is, on a different Hilbert space), 
where x = r is the Euchdean time and = — x is the space variable. Indeed, in that quantization 
scheme, the operator producing Euclidean time translations is e^'^ with H the Hamiltonian 
(generator of time translations). 

This formal analytical continuation has to be made more precise. Consider matrix elements 
of local operators (vac|0(x, t^)|p) with states \p) of P-eigenvalue p. Then, this matrix element 
has singularities as function of p on the positive imaginary axis, and the positions of these 
singularities exactly coincide with the eigenvalues of the Hamiltonian H in the quantization 
scheme where x is Euclidean time. Moreover, the analytical continuation of the matrix element 
towards these singularities gives the matrix element of the same operator in the quantization 
scheme where x is Euclidean time. In relativistic quantum field theory, the singularities are 

branch cuts coming from the measure Yli ^Jpf ~^ ^1 (with p = ^^iPi), and the statement about 
analytical continuation is just crossing symmetry. This claim was also verified to be true in the 
free Dirac theory on the Poincare disk [H], where the singularities are poles and the residues 
must be taken. 

In the case of present interest, our claim is that the "analytical continuation" of the Hilbert 
space C is nothing else than the Hilbert space Tip ol quantization on the circle. This was verified 
explicitly in [2] in the free Majorana theory. 

This claim is made relatively clear by comparing the finite-temperature form factor expan- 
sion (j4.6|) and the expansion in the quantization scheme on the circle (13. 4p . which must agree. 
The analytical continuation we talked about is obtained by shifting the contours of the rapidi- 
ties 9j by the amount ejZ7r/2: then the exponential factors of (14. 6p and of (j3.4p indeed agree, 
under the identification x = t, t = ix. This displacement of the contours can be performed 
while keeping all integrals convergent: we impose x > 0, \x\ > |t| (space-like region) keeping x 
and t fixed, and make only at the end the analytical continuation t = ix. Hence, we see that 
keeping the integrals convergent leads to keeping the operators time-ordered in the quantization 
on the circle. Note that it is here that the condition x > becomes important: the analytical 
conditions that define the finite-temperature form factors will be seen below as consequence of 
this analytical continuation, and depend, at least for twist fields, on our choice of sign of x. 
A different choice of sign would have required a shift in a different direction, and would have 
imposed different conditions on the finite-temperature form factors. 

Of course, the series themselves must agree, but it is natural to assume that they agree 
individually for each term with a fixed number of excitations (at least this can be expected for 
integrable models, where there is no particle production - see section [6] for a discussion). Then, 
these terms will agree for all local fields if the following conditions are satisfied: 

• The factor 

can be written as 

Pai,...,a^,{Ol, ...,9k) F^l...,aki0l^ ■■■,0k; P)F^ak,...-ai{Ok, ■■■,9i;(3), 

where F^^^^ ^^^{9i, . . . ,9k; (3) does not have a pole in the strip lm.{9j) G [0,ej7r/2] for all 
local fields O, and the measure Pai....,a^ (^i, . . . , 9k) has poles at 9j = qj + ej i7r/2 for various 
real qj (more precisely, there are sets = {qi, . . . , qk} for which the measure p has poles 
at such positions, choosing an order: for instance, the pole as function of 9i, whose residue 
has a pole as function of 92, etc.). 
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• We have 

k k ^ 2 ' 

nia^ cosh Qj = Eni,...,nk, ^ "1^^ sinh qj = ^ ^3^' 

where the sets {qj} are the sets sj'^'' for the choice ei = -- - = efc = +) ^-nd are are in one-to- 
one correspondence with the possible configurations of numbers {ni, . . . ,nfc} (there may 

k 

be ambiguities in this correspondence). The values ^ -^S^ must be in Z or in Z -|- ^ (or, 

i=i 

in general, Z + u>), in order to implement the correct quasi-periodicity conditions; this is 
where our finite-temperature form factor expansion (j4.6p is made to agree with the KMS 
identity. 

• The quantities Fai,...,aki^i^ • • • ' ^k', P) for = Qj + ivr/2 and for ei = • • • = e/; = -|-, times 
the square root of the residue of Pai,... ,0^(01, ■ ■ ■ , Qk) at these same values, are proportional 
to the matrix elements ^(vaclOIni, . . . ,nk)i3- 

Note that for the case of one particle, one should recover the energy spectrum En^a = 
ma cosh qn with ma sinh g„ = and nGZif/a = 0,nGZ+iif/a = l. This indicates that 



\ _ — fag— tpma cosh. 6 

and that ¥^{9) = Fa{B). 

Note also that in the free Major ana theory, one simply has 

Pai,...,Ofc(^li ■ ■ ■ ,9k) 



k 

]^ (l _|_ g-ej/3mcosh6lj -j 



with F = F for all excitation numbers, and this indeed reproduces the right energy levels in the 
quantization on the circle as well as the correct matrix elements [2]. In fact, for any free theory 
we have F = F. 

It is now possible to understand the prescription ()4.8p for deforming the integration contours 
in order to avoid possible kinematic poles in the finite-temperature form factors. Indeed, by the 
principles above, the finite-temperature form factor expansion is really an analytical continuation 
of the sum representation of two-point functions coming from the quantization on the circle. 
Hence, it is natural that integration contours be defined in the complex plane to avoid kinematic 
poles, and the direction of the deformation is exactly the one giving the proper correspondence 
between the expansions ()4.6p and ()3.4p . 

It is also possible to understand the restrictions on the directions of the cut of semi-local 
operators, as explained in the second comment after ()4.6p . Indeed, the finite-temperature form 
factor F^^^ ^^{9i, . . . ,9^', P) is, in a sense, the analytical continuation of the matrix element 
(vacjOiIni, . . . ,nk)/^ in the quantization on the circle, which describes "one half" of the two- 
point function. In a path integral formulation, this matrix element corresponds to a path 
integral on the half-cylinder, say, r > 0, x G [0,/3], /3 = with some boundary condition at 
r = (the excited state) and some asymptotic condition at r ^ oo (the vacuum). But since by 
construction the function F^^^^ ^^^{9i, . . . ,9k][3) has no "knowledge" of the other operator O2 of 
the two-point function, it always "stand" in the natural sector given by the trace, even if this 
sector is modified by the cut emanating from O2 in the actual correlation function. Hence, it 
is important that the cut associated to O2 does not change this sector, that is, that it does not 
affect the conditions at r ^ 00 neither those in the region of r present between the position of 
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Oi and that of 02- Similar arguments apply to the function F_^^^ _^_^{9k, ■ ■ ■ ,6i]P), and this 
shows that the cut of Oi must be towards the right, and that of O2, towards the left, when both 
are operators associated to twist fields. 

To be more precise, if the cut associated to O2 does affect the sector in which Oi stands, 
then the only way to provide this information is by modifying the choice of the discrete values 
of rapidities in -^^^,,,^0^.(^1, ■ ■ ■ ,Ok\P) that will form the states on the circle; that is, to modify 
the analytic structure of the measure p. This is indeed what is expected to happen: finite- 
temperature form factors of twist fields should have an analytical structure that provides the 
appropriate shift of the poles of the measure in order to produce the change of sector. Then, 
we actually expect this to be enough information for the finite-temperature form factor of 
Oi to be in this different sector whenever Oi is local with respect to the fundamental fields 
corresponding to the particles involved. In that case, the cut of O2 can indeed be in any 
direction, the expansion ()4.6p will stay valid. Otherwise, if Oi is itself semi-local with respect 
to the fundamental fields, then a change of the analytic structure of the measure is not enough, 
hence the cut of O2 must be in opposite direction. This phenomenon is indeed what is observed 
in the Majorana theory. 

We can circumvent the restriction on directions of the cuts of the twist fields by "twisting" 
the construction; this is done in the next section. 



5 Twisted constructions 
5.1 The twisted inner product 

The construction of the previous section can be fruitfully modified when there is a U{1) inva- 
riance (or sub-group thereof) in the theory, by changing the quasi-periodicity properties of the 
fundamental fermion field in imaginary time. We still consider the space C of endomorphisms 
of 7i, but now as a Hilbert space with a different inner product structure: 

B)^ = Tr(e-/^^+^-Q) " ^'-'^ 

where Q is the Hermitian conserved charge associated to the U{\) symmetry. Then we again 
have 

{A,B)l = {B,A)^. 

Now, we can still consider, in order to have a basis, the set of operators (j4.2p . We have 

where q{a) is the charge^ of the excitation a. Again with diagonal scattering, using this it is 
possible to write down all inner products ()5.ip of the operator (j4.2p in a simple way: 



ii Y — ( — l)fo.i Q2-7riujq{ai) g-ei/3m<i . cosh ' \ ■ J 

i=l ^ ' 

Here, we have assumed the ordering (j4.ip for both members of the inner product. 



^In particular the charge of the fundamental fermion, with e = +, is 1, and in general we have q{a,+) — 
-q(a,-). 
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Note that simple "crossing" relations hold again for operators in EndijL(£): 

^ (.-^ai,...,ak,ce'. 

{9u...,9k,9l),cPL{A)D^,_^, 



,ei_,)) {ei^e.yi) (5.3) 



and 



(6*1, ...,9k), (pR{A)Da'.. 



a'A'^l, ■■■■,^1 



= (^i, ^1, ^fc) , <I^r{A)D^,^_^>{9'^,...,9[)) {9[^9iyi). (5.4) 

5.2 Twisted finite-temperature form factor expansion 

Inverting (j5.2p . we now have the twisted finite-temperature form factor expansion (as before, 
we assume that x > 0): 



{{Oi{x,t)02{m))% 

oo „ 



J- tj {irria a; sinh 6j —irria t cosh ) 



{lm(e,)=e,0+} 



cosh < 



X 



^FZ\...,o.SGi,...,9u;PrF% 



(5.5) 



ife,...,-Ol V^K' • • • ) ' 

where we have defined twisted finite-temperature form factors as the normalised matrix ele- 
ment^: 



= /''.(l_(_l)/.,e2.i..gK)e-./3-..coshe,)J (Ot(o,o),2?,,_,J0i,...,0„))^. (5.6) 

and used the identity 



{^F^^...,adOl, ...,9h(3)y= -F%__^^{9k, ...,9v,P), (5.7) 

which essentially follows from (j6.2p below. Also, we have symmetrised over the orderings of 
rapidities. The quantity AS is non-zero whenever Oi is a twist field a":^^ , and is now given by 



AS 



1 

-+UJ + UJI 



S 



1 

2^^ 



s 



+ UJ — u>i 



iVi = +), 

im = -), 



where if[z/] is the energy of the vacuum |vac^)^ (see the discussion around (j3.9p ). Again, 
when both Oi and O2 are semi-local with respect to the fundamental fermion field, the finite- 
temperature form factor expansion (14. 6p is valid only when the cut of Oi extends towards the 
right (positive x direction) and that of O2 extends towards the left (negative x direction) . This 
is justified in the same way as before, through the relation between (15. 5p and a form factor 
expansion on the circle (j3.4p . Note that with this prescription on the directions of the cuts, one 
now produces the correlation functions 



/3(vac 



i+^+^JOi(x,r)O2(0,0)|vaci^_^^ 



)/3 



with X = T and t = ix. With oJi = —u!2, it is now possible to have the NS vacuum by choosing 

LO = —UJl. 



^ Again, note that the function J'ai,...,ct„(Si, . . . ,9n\(5) has no delta- function contributions at colUding rapidities. 
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6 Formal structure and a generalisation of CFT's mapping 
to the cylinder 

6.1 The space C as a Fock space, and physical interpretation 

We have seen how the space C of operators on TC can be used to obtain infinite series expressions 
of correlation functions (j4.6p and (j5.5p . The important objects are the finite-temperature form 
factors (j4.7p . or the twisted version (j5.6p . which are certain matrix elements on the space C It 
will be convenient, here, to introduce the normalised operators 

k ^ 

dai,...,afe(6'l, . . . , 6'a;) = N Dai,...,ak{Sl, ■ ■ • , ^fc) (6.1) 

with, in the general twisted case, 

9a{9) 



Then, we have 

(^1, ■■■,Ok;P) = {OHo, 0) , . . . , Ok)) 

and the "crossing" relations (j4.4p . (j4.5p change into 

(dai,...,afc(6'l, ■■■,0k) , <AL(^)c^a;,...,a;(^'l> • • • j^';'))^^ 

= • • . • • • (^^ / Vi) (6.2) 

and 

{da„...,a,{ei, ■ ■ ■,9k), <t>R{A)d^'^_^,{9[, ■ ■ ■,9[))^ = e2W";)e-^'i/^'"«; -^^^'i 

X (da;,ai,...,a,(^'i, 9i,..., 9k) , </'R(^)da'„...,a;(e2, ■ ■ ■ , 9[))^ {9[ / 0, Vi). (6.3) 

In order to describe in a convenient way the space spanned by dai,...,^^. (^i, . . . ,9k), we intro- 
duce the following operators acting on C: 

4(0),Z,(0) GEnd(£) 
with the following properties: 

Zi{9)l-H = da{9), ZU^l)da,,...,ad^2, ...,9k) = (^1 , . . . , 9k), Z^{9)ln = 

and satisfying the following exchange relations: 

Za■^{9l)Za^{92) — Sai,a2{9l — ^'2)^^02 (^2)^01 (6*1) = 0, 

ZU9^)ZI{92) - 5«„«,(0i - 92)Zl^{92)Zi^{9i) = 0, (6.4) 
Z^,{9i)Zl{92) - 5„„„,(02 - 9i)Zi^{92)Z^,{9i) = -J—6^^^,^6{9, - 92), 

9ai[9l) 

where (Jq^^oj = ^ai,a2^ei,e2 and where 

e {n\ _ j ^a.i,a2{9) (ei = £2), 

y Sa2,aA-9) (ei = -£2). 
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The space 7i is then seen as a Fock space over the algebra ()6.4p . with vacuum vector 1-^ an- 
nihilated by Za{9). The algebra (j6.4|) has exactly the structure of the Zamolodchikov-Faddeev 
algebra (j2.ip with diagonal scattering, but with twice as many particles. The physical interpre- 
tation is that the "states" c?«i,...,a^ (^i, • • • ,6^) correspond to configurations of stable "additional 
particles" (e = +) and "missing particles" or holes (e = — ) in a thermal bath (we will call 
both "excitations"), both created by the operators Za{0), and scattering through the matrix 
'S'ai, 02(^1 ~ ^2)- They are stable, since the states with n excitations have no overlap with those 
with n' ^ n excitations. It is both this stability and the fact that matrix elements should 
have nice analytical properties that suggests that the correspondence between expansions in the 
quantization on the circle and finite-temperature form factor expansions holds individually for 
every term with a given excitation number. Only in integrable quantum field theory can we 
expect these two properties together. 

In order to have a better picture of the "particle" and "hole" states that we are discussing, 
one should recall that we defined the in-states of 7i by multiple action of Za{9) on |vac) with 
ordered rapidities, decreasing from left to right. When this order is not satisfied, one has an 
"intermediate" state, which is a useful concept only in integrable quantum field theory. It 
corresponds to having wave packets ordered such that some interact in the far past (like for 
out states), while others interact in the far future (like for in states). Essentially, the order of 
the operators Za{0) acting on |vac) corresponds to the order of the wave packets themselves 
along the x axis, when time is taken to go "upward" (in the positive y direction). When we put 
an operator Za{0) inside the finite-temperature trace (evaluating the trace in the in basis, for 
instance), we take situations with various numbers of in particles at various ordered rapidities, 
and put an additional wave packet far to the left, generically producing "intermediate" states. 
This is the sense in which the operator Za{0) corresponds to an additional particle in a thermal 
bath. Similarly, Za{0) is taking away a particle from the thermal bath, by first bringing its 
wave packet far to the left. Recall that it is because we keep the order of the wave packets fixed 
while varying the rapidities that matrix elements of operators on Ti in "intermediate" states are 
meromorphic functions. With the previous discussion, this lends support to the fact that the 
basis -Dai,...,afe(^i) • • • )^fc) of C should produce matrix elements with nice analytical properties. 

It is also interesting to note a nice physical interpretation for the expected kinematical poles 
in the finite-temperature form factors, occurring at colliding rapidities when they are associated 
to opposite signs of e. These poles can be seen, in the finite-temperature form factor expansion, 
to lead to powers of the time variable t, instead of exponential factors. This corresponds to the 
fact that a particle and a hole can annihilate and re-form themselves at arbitrary large time 
differences without cost in energy. For a given excitation number, various powers in t will occur, 
which can more or less be made in correspondence with various simultaneous annihilating and 
re-creating of particle-hole pairs. It it these processes that make the computation of large-time 
dynamics from finite-temperature form factors difficult: a re-summation of these powers of t is 
necessary. 

The fact that the Hilbert space C is the same as an ordinary Hilbert space with twice as many 
"particles" does not mean that the finite-temperature theory is the same as a zero-temperature 
one with such particles. Indeed, another ingredient that defines a theory is the set of local opera- 
tors (in particular, the energy density), and this set looks very different on a finite-temperature 
Hilbert space. In fact, it would be very interesting to study the structure of the energy density 
on C. 

6.2 A tool for evaluating traces and mapping to the cylinder 

This is again a generalisation of a concept introduced in [2]. It is possible to construct the left 
action of Z^ through the operators and Z^ (the right action is more complicated, from the 
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definition of the operators Za and Za, but it is sufficient to consider the left action). Indeed, 
one can verify that the mapping 

MZaiO)) = Zi{9)g^{6) + Z^{0)g^{0) (a = (a, -e) for a = (a, e)) (6.5) 

is an isomorphism of the algebra satisfied by Za[9): 

using the property 502,01(^2 — ^'i) = Sai,a2{^i — ^2) as well as ga{0) + (— l)^~'^"5a(^) = e^~^°-. 
From this one can conclude the equality 

{inAL{Zo.M) ■ ■ ■ </>L(^,,(0fc))iH) = Hh^z^AOi) ■ ■ ■ Zo^M)) 

by bringing on both sides, for instance, all factors with e = + to the right, by evaluating the 
resulting right-hand side explicitly using the inner product (15. 2p . and by evaluating the left- 
hand side using the representation (|6.5p . Then, restrict all rapidities numbered j + 1, . . . , /c to 
be different from one another, and restrict all rapidities numbered 1, ... ,j — 1 to be different 
from one another. The left-hand side evaluates to 

{Zaj^i{Oj-l) ■ ■ ■ Zai{9i),4>L{Zaj{0j))Za^^-i^{6jj^i) ■ ■ ■ Za^{Ok)), 

whereas the right-hand side evaluates to 

{Zaj^A^j-l) ■ ■ ■ Zai{9i), Zaj{ej)Zaj^^{9j+i) ■ ■ ■ Za^{9k))- 

The equality of these two expressions shows the equality for all matrix elements, hence 
shows ()6.5p . 

Prom relation (j6.5p . one can write explicitly products Za^[9i) ■ ■ ■ Za^{6k) as linear combi- 
nations of the basis elements iia'^,...,Qj(^i, • • • ,9[) (for / < k and where primed variables form 
a subset of un-primed variables), using 

Za., {9l)--- Za,M =MZaA(^l))--- MZa, {0k))lH- 

In particular, using cyclic properties of the trace, one can relax the restriction of having different 
rapidities in ()6.2p . and one obtains 



= {da[,...,a'i_i(^l^ ■ ■ ■ ,Si-i),4>l{A) daj^^...,ak,a[{Gl, ■ ■ ■ ,0k,0[)) 

v-U ^a,a'^{0j~0'l) 
+ 2^ •S'oj.aj+i {Oj - 9j+i) ■ ■ ■ Sa^^ak {(^j,&k)^^ 

X (4;,. (^'1, • • • , Cl), 4>l{A) ...,9j,..., 9^)), (6.6) 

where the hat means that the variable is missing. 

It is now possible to generalise, in some sense, the concept of "mapping to the cylinder" 
that can be used in conformal field theory in order to evaluate finite-temperature correlation 
functions. In conformal field theory, one has a mapping (of vertex operator algebras) O ^ O 
such that correlation functions of O on the cylinder are equal to correlation functions of O on 
the plane (for instance, for the stress energy tensor one has T i— > -|- z~'^ where w is 

the coordinate on the cylinder, z = e^^"'/^ is the coordinate on the plane and c is the central 
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charge). We can generalise this here at the level of form factors. We seek a linear map Q from C 
to C such that 

(lH,<AL(f^^)c^(ai,+),...,(afe,+),(a; ,-),...,« • • ■,0k,0l, . . .,9[)) 

= a'„...,a',i(^'l, . . . , el\A\ei, ek)a,,...,a, {0^ + 0, ^ij), (6.7) 

where A £ C. This requirement is inspired by the fact that 

^11^(1^1 4>L{A)d(^ai,+),...,iak,+),{a[~),...,{a[~){01, ■ ■ ■ , ^k, O'l, . . . , O'l)) 
= ai,...,a;(^l, • • • , e'lim, ek)a,,...,a, (0- + 0j Vi, j) 

thanks to (j6.2p . Note that this in fact completely fixes the map (j)L°^, thanks to (j6.6p . and by 
injection the map Q, if it exists. 

In order to describe the map 0,, it is convenient to consider elements of C that have simple 
expectation values on 7i. We consider products where operators Za{9) are normal-ordered with 
respect to the vacuum in 7i: operators Z"'{9) are placed to the right of all operators Za{9), 
without taking any delta- function term (but taking all S'-matrices involved in the exchanges). 
We will denote the normal-ordering of A by the standard '.A:. The set of all normal-ordered 
operators spans C. On the other hand, there is a natural normal-ordering with respect to the 
vacuum l-j^ in C: it is the one whereby operators Z are placed to the right of operators . We 
will denote the normal-ordering of ^ G End(£) by *^*. Then, it is simple to see that for all 
normal-ordered A (that is, :A: = A), 

MA): (6.8) 

can be put in place of (pLi^iA)) in (j6.7p in order to have equality. Indeed if A contains exactly 
I operators of type Za and k operators of type Z^, then the equality is clear; otherwise, both 
sides are zero, hence the equality still holds. However, this does not yet show that the map 0, 
exists. 

Now, let us choose a basis in C, with normal-ordered elements Ai composed of products of 
finitely many operators Z^. Certainly, the set (j)L{Ai) is not a basis in End(£). However, if we 
are to project from the left with 1^^, then we do obtain a basis: {l-n, (j)L{Ai)-) is a basis in the 
dual of C. Similarly, the space of operators t</'L(^j)* gives another basis in the dual of C when 
projected from the left with !■}{. Hence, we have the change-of-basis relation 

(In, ■)=Y1 cpUAj) ■), (6.9) 

j=0 

where the sum is finite. Then we find, from the discussion around (j6.8p and recalling that 
■Ai'. = Ai, 

QAi = y^^MjjAj. 

j=0 

The change of basis (j6.9p could be calculated explicitly, but there is a nice way of expressing 
it. Suppose we can find an operator T in End(£) such that 

Tin = 0, [T, Z^ie)] = 0, [T, Ziie)]gUO) = ^^(0)g(,,„)(0). 

Then one can verify that e^J</>L(:^:)*e~^l7^ = :A:, hence that 
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But using the algebra (j6.4p . one can see that 

T = / deY,2^a,-){0)Zia,+)i0)9ia,-)iO) (6.10) 

has the right properties. Hence we have found an exphcit expression for the map, acting on the 
space of operators on 7i, that transform finite-temperature form factors into zero-temperature 
form factors. This is the generahsation of the concept of "mapping to the cyhnder" in conformal 
field theory. Finite-temperature form factors can then be calculated using 

(1h,<Al(^) d(ai,+),...,{ak,+),{a[ -),... ,{a[-){Gl, ■ ■ ■ ,Ok,0'u ■ ■ • i^'l)) 
= a;,...,a;(^'l) • • • (e^^) l^li • • • ,(^k)ai,...,ak- 

It is important to note, however, that the operator does not act on the initial Hilbert space Ti, 
but rather on the space of operators acting on it C In the quantization on the line, this is not 
isomorphic to TC. In conformal field theory, one usually thinks about the quantization on the 
circle around a fixed point (radial quantization), and by the operator-state correspondence, one 
then has an operator acting on the Hilbert space that performs the mapping to the cylinder. 

The action of on ^ S £ can be made more explicit using ()6.10p . Indeed, taking A = 
Zai (^i) • • • Za^. {9k), we can write 

e^iZ^M • • • ZaM)) = e^MZa^iOi)) • • • <PL{ZaM))'^H 

= e^MZa^ {Oi))e-^ ■ ■ ■ e^MZaM)e-^lH 

and use 



e 



4>L{Za{9))e-' = MZa{9))+2a{9)gia,-){0), 



Za{9)MZa'{e')) = Sa,a'{0 " 0')0L (^a' (^') (^) + Sa,a'6{e - d') 

to bring all the 2^(9) to the right, annihilating 1-^. This leads to a natural generalisation of 
Wick's theorem, whereby A is written as j4 + operators where more and more contractions 
have been performed, the contractions being given by 

contraction of Za{6) with Za'{0') = 5(a,-)(^)'^a,a'(^ ~ ^0- 

If the S'-matrix Sa,a'{G ~ G') is equal to ±1 (free models), then this gives the standard Wick's 
theorem, and in particular it can be applied to linear combinations of the type j d9Za{9)fg 
as well. In free models certain linear combinations of this type indeed represent local fields, 
and this immediately leads to the phenomenon of "mixing" that was described in [2] (using 
slightly different arguments). Of course, local fields in interacting models and twist fields in 
general are not simply such linear combinations, but rather are sums of operators with more 
and more factors of Za{9) (since many-particle form factors are non-zero). Hence all term will 
contribute to any given finite-temperature form factor, and it is a non-trivial matter to re-sum 
these contributions. 

7 Results in the free massive Majorana theory 

7.1 Free massive Majorana fermions 

The free massive Majorana theory with mass m can be described by the action 
A = i I (fx{—ip{dx + dt)^ + ip{dx — dt)i) — mipip). 
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It is a model with only one particle, and with only Z2 internal symmetry, described by a change of 
sign of the fermion fields. In particular, the fields ijj and iJj are both real (hence the corresponding 
operators in any quantization scheme are Hermitian). The quantization on the line is simple to 
describe. Fermion operators are given by: 



where the mode operators a{9) and their Hermitian conjugate a\9) satisfy the canonical anti- 
commutation relations 

{a\9),a{9')} = 6{9-e') (7.1) 
(other anti-commutators vanishing) and where 

pg = msinh^ , Eg = mcosh9. 
The fermion operators satisfy the equations of motion 

— 1 TTl — 

di){x, t) = - (9x + dt)tp = 

— 1 — 771 

d^{x,t)^-id^-dt)^P = -^ (7.2) 

and obey the equal-time anti-commutation relations 

{ip{x,t),^p{x',t)} = 6{x-x'), {iP{x,t),^P{x',t)} = 6{x-x'). (7.3) 

The Hilbert space TC is simply the Fock space over the algebra (j7.ip with vacuum vector |vac) 
defined by a(0)|vac) = 0. Vectors in 7i will be denoted by 

\9i,...,9k) = a\9i)---a\9k)\vac). 

A basis is formed by taking, for instance, 9i > ■ ■ ■ > 9k- This is exactly the construction 
described in Section [21 with only one particle and S{9) = —1. The Hamiltonian is given by 



/oo 
d9m cosh 9 a\9)a{9) 
-00 



and has the property of being bounded from below on TC and of generating time translations: 

[H, ^{x, t)] = -i-g^Hx, t), [H,i>{x, t)] = -igi^x, t). (7.4) 

In the discussions of the previous sections, we also considered quantization on the circle of 
circumference /?. It will be convenient to have the description of this quantization for the present 
model, with anti-periodic (NS) conditions on the fermion fields. The fermion operators evolved 
in Euclidean time r are: 



nG 



an/2 



V^(x,r) = — y= I One— -a;e 

V2L vcosha„ V 
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where the discrete mode operators and their Hermitian conjugate alt satisfy the canonical 
anti-commutation relations 

{al,an'} = Sn,n' (7.5) 

(other anti-commutators vanishing) and where 

27rn / 1\ , , 

Pn = msinha„ = — — n G Z + - , ('-6) 



L V 2^ 

En = mcosha^. 

The fermion operators satisfy the equations of motion (j7.2p as well as the equal-time anti- 
commutation relations (17. 3p (with the replacement ip and ip ip); the latter is simple to 
derive from the representation 



5{x) = Y Yl ^'^""^ 



of the delta-function, valid on the space of antiperiodic functions on an interval of length /3. The 
Hilbert space 7^/3 is simply the Fock space over the algebra (j7.5p with vacuum vector |vacNs)/3 
defined by an|vacNs)/3 = 0. Vectors in 7ii3 will be denoted by 

|ni,. . .,nk)i3 = ■ • • aJ,JvacNs)/3• 
A basis is formed by taking, for instance, ni > • • • > n^. The Hamiltonian (with vacuum energy) 
is given by 

Hf3 = £ns + ^ m cosh On a^an 

and has the property of being bounded from below on TC/^ and of generating time translations: 

d - ^ d ^ 

[Hp, ^(x, r)] = ^V'(x, r), [Hp, t)] = ^V'(x, r) . 

Our discussion was with the NS sector in mind, but it is not hard to perform the quantization 
in the R sector. What will be important for us are relative energies of the NS and R vacua: 

/OO Jfl , , 

— cosh^lnfl + e-"^^™^*^^) , 
-OO 27r V / 

/OO jn , . 

— cosh01n(l-e-"^^'=°^^n , (7.7) 
-OO 27r V / 



where we used the notation of the discussion around (j3.9p . Here, the vacuum energies of both 
sectors were calculated in the same regularisation scheme and e contain terms that are common 
to both. 

It is worth noting that the normalisation that we took is slightly different from the more 
standard normalisation in conformal field theory, that makes the fields and ip not real, but with 
definite phase. With our normalisation, the leading terms of the operator product expansions 
(OPE's) ^/;(x,t)V'(0,0) and iij{x,t)il;{0,0) are given by 

^(x, t)V(0, 0) ~ ' ^(x, t)i;{0, 0) ~ - \ (7.8) 

27r{x — t) 27r(x + t) 



Finite- Temperature Form Factors: a Review 



25 



7.2 Twist fields 

Two fields are of particular importance: they are two primary twist fields associated to the Z2 



symmetry, which we will denote by a and fi as is customarjo, the first one being bosonic, the 
second fermionic. In the sense of quantum chains, the first one is an "order" field, with non- 
zero vacuum expectation value, the second is a "disorder" field, with zero vacuum expectation 
value. As we explained in sub-section 13.31 to each of these fields there are two operators on 7i, 
which makes four operators: a± and //±. They are fully characterised by the leading terms in 
the (equal-time) OPE's that are displayed in Appendix |Al These leading terms are fixed by 
the general requirements ()3.5p and (13. 6p , by our choice of branch which says that when fermion 
operators are placed before the twist-field operators, they are on the same branch no matter the 
direction of the cut, and by the general "field" product expansion that holds inside correlation 
functions: 

^(x, t - «0+)a(0, t) ~ ' u(0, t), ^{x, t - «0+)^(0, t) - ^ a{0, t), 
2\/ttx + iO+ 2\/ttx + zO+ 

ij{x, t - iO+)a(0, t) ~ - , ' .^^ /i(0, t), ij{x, t - iO+)/x(0, t) ~ ^ a(0, t) 

2\/ttx — iO~^ 2vvrx — zO+ 

with branch cuts on x < 0. 

It is worth nothing that the relations of Appendix [X] are in agreement with the Hermiticity 
relations a^. = a± and = ±//-|-. 



7.3 Riemann— Hilbert problem for twisted 

and untwisted finite-temperature form factors 

7.3.1 Untwisted case 

In [2], the (untwisted) finite-temperature form factors (|4.7p of the twist-field operators above 
were shown to solve a Riemann-Hilbert problem of the type found at zero temperature, but 
with important modifications. We repeat here the results. 
Consider the function 

frjiOi, ... ,6k) = F^^^ _^(6'i, ...,6k](i) 

where 0^ is the operator with branch cut on its right (r/ = -|-) or on its left {r] = —) representing 
any twist field: this can be the order field a± or the disorder field fi±, or any of their conformal 
descendants (that is, fields which reproduce conformal descendants in the massless limit). Con- 
formal descendants include space derivatives, as well as other fields related to action of higher 
conformal Virasoro modes on twist fields. A way of describing such descendants is by taking 
the limit x ^ of the finite part of the OPE O{x)a±{0) or O{x)fi±{0), where O is any bosonic 
operator formed out of normal-ordered products of fermion operators. 
The function / solves the following Riemann-Hilbert problem: 

1. Statistics of free particles: / acquires a sign under exchange of any two of the rapidity 
variables; 

2. Quasi-periodicity: 

^(01,.. .,9j + 2m, ...,9k) = -f{Oi, ...,0j,.. .,9k), j = l,...,k; 

*In the present section, the symbol a does not denote a generic twist field, but rather the primary twist field 
as described here. 
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3. Analytic structure: / is analytic as function of all of its variables 6j, j = 1,... ,k ev- 
erywhere on the complex plane except at simple poles. In the region lm{6j) G [—in, in], 
j = I, . . . ,k, its analytic structure is specified as follows: 

(a) Thermal poles and zeroes: frjiOi, . . . ,6^) has poles at 

m 

Oj=an-'n—, n£Z, j = l,...,k 
and has zeroes at 

ZTT 1 
0j = an - 7]—, neZ+-, j = l,...,k, 

where a„ are defined in (j7.6|) (and, of course, we use this definition for any n); 

(b) Kinematical poles: f^{6i, . . . ,6^) has poles, as a function of 6^, at 9j ± in, j = 
1, . . . ,k — 1 with residues given by 

f^^^'^ • • • ' ~ n i_e-/^^«. e^^e;^^ • 



In order to have other finite-temperature form factors than those with all positive charges, 
one more relation needs to be used. We have: 



4. Crossing symmetry: 

F?,:..,,,,...,e,{ei, ...,9,+in,...,ek;l3) = iFf^l,^_,._,^{9i, 



The name "crossing symmetry" is inspired by the zero-temperature case (and it is not to be 
confused with the simpler "crossing relations" introduced in (j4.4p . (j4.5p . (j6.2p . (j6.3p ). To make 
it more obvious, define the functions 

fr,{e[, e[\e^, ...,9^) = id+_+{e[, 9l), <Al(o,(o, o)) ^^+,...,+(^1, . . . , 9k). (7.9) 

These are in fact distributions, and can be decomposed in terms supported at separated rapidities 
9[ 7^ 9j, yi,j, and terms supported at colliding rapidities, 0- = 9j for some i and j. We will 
denote the former by f^''^-{9[, 9'i\9i, 9k), and the latter by /^°"-(^^ • • • , &k)- 

Under integration over rapidity variables, the former gives principal value integrals. Recalling 
the property (j6.2p . we have 

/^^P'(6'i, . . . ,9'i\9i, . . . ,9k) = -^+,^..,+ -(^1, • • • ,9k, 9'i, ...,9'i;/3) 

for {0'^ / 9j Vi G {!,...,/}, j £ {1, . . . , k}), where on the right-hand side, there are k positive 
charges (+), and / negative charges (— ). Analytically extending from its support the distribution 
/^°^' to a function of complex rapidity variables, crossing symmetry can then be written 

/^^P-(e;, . . . , 9[\9u ...,9k + in) = i/^<'P-(^l, . . . , 9'i, 9k\9i, 9k^i), 
/-P-(ei, ...,9[ + in\9i, ...,9k) = ifr,''''-{9[, 9'i_,\9i, 9k, 9[), 

which justifies its name. 

It is worth mentioning that the distributive terms corresponding to colliding rapidities satisfy 
a set of recursive equations: 
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5. Colliding part of matrix elements: 

Jri [f^lj ■ ■ ■ i^lWli ■ ■ ■ i^k) 

Ik f3Eg. 

i=i j=i 1 - e ' 

Note that the colliding part vanishes in the limit of zero temperature, (3 — > oo. Finally, it 
is instructive to re-write the distribution fr^{9'i, . . . ,9'i\9i, . . . , 9^) as an analytical function with 
slightly shifted rapidities, plus a distribution, using the relations 

^ ±z^5(^)+pfiV (7.10) 



9^m 



where P means that we must take the principal value integral under integration. Defining the 
disconnected part . . . ,9[\9i, . . . ,9^) of the matrix element (j7.9p as 

fr,{9[,.. .,9[\9^,...,9k) = fn'^^'iO'i - ??iO+, ...,9'i- 7?iO+ I 01, . . .,9k) 

+ f^~{9[,...,9[\9i,...,9k), 

where again we analytically extend from its support the distribution f^^^' to a function of 
complex rapidity variables, we find that the disconnected part satisfies the recursion relations 

f disconn. rnf Q> I Q Q \ 

Jn [f^l, . . . . . . ,^k) 

I k 

= E (1 + ^"'''''^ ) ^(^^ - ■■■Jl---,0l\9i,...J„..., 9k). 

i=i j=i 

Note that the factor (1 + e ) 5{9[ - 9j) appearing inside the double sum is just the overlap 
{d+{9^),d+{9j)), so that the equation above can be naturally represented as a "sum of discon- 
nected diagrams." This equation is, in fact, consequence of the general relation (16. 6p . 



7.3.2 Twisted case 

The twisted case was not considered in [2], but can be obtained from the same arguments. 

There is no U{1) invariance, but we can still twist by the Z2 symmetry. Hence, we consider 
twisted finite-temperature form factors ()5.6p with uj = The derivation of |2j for the Riemann- 
Hilbert problem can easily be adapted to this case, and the results are as follows. 

Consider the function 

frj{9i, . . . ,9k) = 2F^J ^^(^1, . . . , 9k; P), 

where is the operator with branch cut on its right (r/ = +) or on its left (r/ = — ) representing 
a twist field. The function / solves the following Riemann-Hilbert problem: 

1. Statistics of free particles: / acquires a sign under exchange of any two of the rapidity 
variables; 

2. Quasi-periodicity: 

^(01,.. .,9j + 2i7T, ...,9k) = -UOi,.. .,9k), j = l,...,k; 

3. Analytic structure: / is analytic as function of all of its variables 9j, j = 1, . . . ,k every- 
where on the complex plane except at simple poles. In the region lm{9j) £ [—in, in], 
j = 1, . . . , k, its analytic structure is specified as follows: 
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(a) Thermal poles and zeroes: fri{Oi, . . . ,6^) has poles at 

ivr 1 
Oj =an-r]—, neZ+-, j = l,...,k 

and has zeroes at 

6j = an-r]—, n G Z, j = l,...,A;; 

(b) Kinematical poles: /^(^i, . . . , 0^) has poles, as a function of 9^, at 6j it ivr, j = 
1, . . . , /c — 1 with residues given by 

Again, in order to have other finite-temperature form factors than those with all positive 
charges, one more relation needs to be used. We have: 

4. Crossing symmetry: 

-^FF,:..,e„...,e,{Ou ...,9j+i7r,...,9k;P)=i ^Fe?:..,-.,,...,.. (^1, • • • , ^i, • • • , Ok-, 13). 
Moreover, matrix elements 

f^{9[, 9l\9i, ...,9k) = . . . , 91), MOviO, 0)) d+_+{9i, ...,9k). 



can again be decomposed in terms supported at separated rapidities 9[ ^ 9j, (which 
give principal value integrals under integration), and terms supported at colliding rapidities, 
~ % fo'^ some i and j, denoted respectively by f^'^^'{9'i, . . . ,9i\9i, . . . , 9^) and /^"^^■{9[, . . . ,9'i\ 
9i, ... ,9k). Recalling the property (j6.2p . we have 

f^'^-{9[, . . . ,9[\9i, . . . ,9k) = 2F^^^ _,_ . . . ,9k,9[, . . . , 9[;f3) 

for (6*^ 7^ 9j Mi G {\, . . . ,1}, j G {1, . . . , k}), where on the right-hand side, there are k positive 
charges (+), and / negative charges (— ). The distributive terms corresponding to colliding 
rapidities satisfy again a set of recursive equations, now modified by the twisting: 

5. Colliding part of matrix elements: 

j^coU. ( nl qI \ Q Q \ 

Jri \Pl, ■ ■ ■ i^l]^!-, ■ ■ ■ ,^k) 

= E '^(^^ - ...A,---A\0i....,9„..., 9k). 

i=i j=i 1 + e J 

Finally, we can again re-write the distribution f^{9[, . . . ,9i\9i, . . . ,9k) as an analytical func- 
tion with slightly shifted rapidities, plus a distribution, using the relations (|7.10p . Defining the 
disconnected part /^'^™'^''-(6li, . . . ,9l\9i, . . . ,9k) oi the matrix element ([72]) as 

f^i9[, ...,9'i\9u...,9k) = /^'^P-C^l - ??iO+, ...,9l- 7?iO+ \9i,...,9k) 

I f disconn. (nf Q> \ Q Q \ 

+ Jri [t^l, ■ ■ ■ ■ ■ ■ ,^k), 

where again we analytically extend from its support the distribution f^^^' to a function of 
complex rapidity variables, we find that the disconnected part satisfies the recursion relations 

I k 

= E E(-l)''"'"'"' (1 - ^"'''''^ ) ^(^i - ■■■,0i,---,Ol\9i,... ,§„..., 9k). 

1=1 j=i 
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7.3.3 Other local fields 

It is worth noting that points 1, 2 and 4 are in fact also valid for fields that are local with 
respect to ip and ip. The analytic structure, point 3, for such fields, is much simpler: the finite- 
temperature form factors are entire functions of all rapidities. In fact, the finite-temperature 
form factors of tp and ip themselves are exactly equal to their zero-temperature form factors, and 
for other fields, a phenomenon of mixing occurs, as described in [2] and as can be calculated 
using the techniques of Section [6l 

7.3.4 Differences with zero-temperature Riemann Hilbert problems, 
and some explanations 

There are three main differences between the Riemann-Hilbert problems stated in this sub- 
section, and the Riemann-Hilbert problems solved by zero-temperature form factors, reviewed 
in Section [2j First, there are, in the former, so-called "thermal" poles and zeroes. They are in 
fact consequences of the semi-locality of the operators with respect to the fundamental fermion 
fields, and play the role of "changing the sector" of the excited states when integrals are deformed 
to reproduce the form factor expansion in the quantization on the circle. Indeed, they displace 
the poles of the measure in order to reproduce the right set of discrete momenta. 

Second, the kinematical residue has an additional factor. This factor, in fact, is closely related 
to the presence of the thermal poles and zeroes. 

Finally, there is a subtle but important difference: the quasi-periodicity equation has a sign 
difference. Essentially, the quasi-periodicity equation that we have at finite-temperature is 
exactly the one we would have at zero-temperature with fields that are local with respect to 
the fermion field. This difference is again due to the presence of the thermal poles and zeroes. 
More precisely, in the limit of zero temperature, the finite-temperature form factors converge 
to the zero temperature one only in the strip Im(^) g] — 7r/2,7r/2[. At the lines lm(0) = ib7r/2 
(the sign depending on which excitation type and form factor we are looking at), there is 
an accumulation of poles and zeroes that produces a cut. The quasi-periodicity equation of 
zero temperature comes from the analytical continuation through this cut. Note that it is this 
analytical continuation that recovers rotation invariance in Euclidean plane, an invariance which 
is broken by the cylindrical geometry at finite temperature. 

We would like to mention, in relation to the breaking of Euclidean rotation invariance, that 
yet crossing symmetry, point 4, is valid. It is in fact a consequence of the fact that the defor- 
mation of the contours, as explained in Subsection 14.41 should give residues at the poles of the 
measure occurring in ()4.6p . These residues come from two contributions: the contribution of the 
displaced 9 contour associated to e = -|-, and that associated to e = — . That these two contribu- 
tions should give a residue impose certain conditions on the value of the finite-temperature form 
factors: they should correspond to contours in opposite direction and on opposite sides of the 
same pole. From this and from knowing that all finite-temperature form factors of the fermion 
fields satisfy crossing symmetry, one concludes that crossing symmetry holds for all local fields. 

7.4 Results for twisted and untwisted finite-temperature form factors 

Again, we repeat here the results of [2], and generalise them to the twisted case. Note that the 
method of computing one-particle finite-temperature form factors by solving the Riemann- 
Hilbert problem with this asymptotic is very similar to the method used by Fonseca and 
Zamolodchikov |21j for calculating form factors on the circle. 

For the order and disorder operators, a± and fj,± the solutions to the Riemann-Hilbert 
problems above are completely fixed (up to a normalization) by the asymptotic ~ 0(1) at 
9 zboo, imposed by the fact that they are primary fields of spin 0. 
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For the one-particle finite-temperature form factor of the disorder operator with a branch 
cut on its right, the solution is 



C{(3) exp 



1 



ooTiO+ 2vri sinh(6l - 6') 



In 



l + e 



-l3Eg 



(7.11) 



for some real constant C{[3). This is in agreement with the Hermiticity of /U+, which gives 
{F^+{9;f3))* = F^+{e;P) for 9 real. Using 



1 



sinh(6' - {6' ± iO+)) 
this can also be written 



±i7rS{e - 61') + P 



1 



sinh(6l - 9') 



^^+{9;P)=C{P)e 



±2 




-PEg 



-l3Ee 



exp 



d9' 



1 



2-Ki- Vsinh(6l - 9' 



In 



l + e 



-PE, 



That this is a solution can be checked by verifying the asymptotic F^{9;I3) 



~ e 



± ■ 



C{I3) as 



cxD, and by verifying that the functions F'^{9\(5) have poles and zeros at the proper 



positions. Positions of poles and zeros are the values of 9 such that when analytically continued 
from real values, a pole at smh{9 — 9') = in the integrand of (|7.1ip and one of the logarithmic 
branch points pinch the 9' contour of integration. The fact that these positions correspond to 
poles and zeros can be deduced most easily from the functional relation 



1 I -f3Eg 

F£+ {9; (3)F^^ {9 ± zvr; P) = ±iC{Pf-^^-^. 



(7.12) 



Note that this implies the quasi-periodicity property 
+ 2m-p) = -F^+{9-(3). 



It is also easy to see that the crossing symmetry relation is satisfied. 

For the operator with a branch cut on its left, one can check similarly that the function 



-iF^*{9;P) 



■{9-p) = F^+{9-m-(3)- 
solves the Riemann-Hilbert problem of Paragraph 17. 3. II with r] 

,OO±i0+ ^Q, ^ 



F^- 



'-C{(3) exp 



-j±iQ+ 2-Ki sinh(0 — 9') 



In 



. Explicitly. 

1 + e-^^' 



1 



(7.13) 



{9; P), which is in agreement with the anti- 



In particular, we observe that (-F±" {9; /?))* = —F^ 
Hermiticity of the operator Note that we chose the same constant C{P) as a normalization 
for both Ft^~ and Fl^'^ . This is not a consequence of the Riemann-Hilbert problem, but can 
be checked by explicitly calculating the normalisation. The normalisation was calculated in [2], 
and is given by 



CiP) 



(7.14) 



where the average ((c))/? was calculated in [28] (the average at zero-temperature (that is, P — > cxd) 
can be found in [55]) and is given by 



1 J_ -1.3 

m 8 2 12 e exp 



{mPf 



d9id92 



sinh 9i sinh 92 



(27r)2 sinh(m/? cosh ^i) sinh(?7T,/3 cosh ^2) 
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X In 



coth — 



where A is Glaisher's constant. Essentially, this normalisation is evaluated by computing the lea- 
ding of {{ip{x,0)fi^{0,0)))f3 as X — > 0+, and the leading of {{fi^{0,0)^{x,0)))f3 as x — > 0~, using 
the form factor expansions; in both cases, it is important to approach the point x = from a 
region that is away from the cut. 

Multi-particle finite-temperature form factors can be easily constructed from the well-known 
zero-temperature form factors (first calculated in [9]), by adjoining "leg factors", which are just 
normalized one-particle finite-temperature form factors: 

F,°-,,(^,,...,^,;/3)=Jt]((a))Jn^t^ JJ tanh , 

where 0+ is 0"+ if k is even, and ^+ if k is odd. The symbol [k/2] equals the greatest integer 
smaller than or equal to k/2. This satisfies the condition on thermal poles and zeroes simply from 
the properties of the leg factors, and it can be verified that this satisfies the quasi-periodicity 
condition and the kinematical pole condition. Point 2 and Point 3b of Subsection 17. 3. H respec- 
tively. Using crossing symmetry, Point 4, it is a simple matter to obtain the formula for other 
values of the charges: 

(n^) J J-K^^ 

Similarly, we have 

where O- is a- if k is even, and fi- if k is odd. 

Finally, twisted one-particle finite-temperature form factors can easily be obtained by solving 
the Riemann-Hilbert problem of Paragraph 17.3.21 as follows: 



Fr{o) 



These functions have the correct analytic structure, they satisfy the crossing symmetry relation 
(point 4), and their normalisation is the correct one that can be deduced from the fact that the 
leading of ('0(x, 0), /x_ (0, 0)) i as x ^ 0"*", and the leading of (/i+(0, 0), '(/'(x, 0)) i as x ^ 0~, are 

2 2 

the same as in the untwisted case. Twisted multi-particle form factors can also be obtained in 
a simple way: 



where 0± is a± if k is even, and /i± if k is odd. 
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7.5 Form factors on the circle from finite-temperature form factors 

As explained in Subsection 14.41 there is a relation between finite-temperature form factors and 
form factors in the quantization on the circle. In the present case of the Majorana theory, this 
relation was written explicitly in [2], and was proven by independent means. A slight extension 
to the twisted case gives it as follows: 



^(ni, . . . ,ni|O(0,0)|ni, . . .,nk)p 

{-) 

\mL J 




X ( + y,...,a„, + y,a^,, + y,...,an, + y;/?J , (7.18) 

where there are k positive charges and / negative charges in the indices of '^F'^, and where an 
are defined in (j7.6p . Here, s is the spin of O. This formula is valid for any excited states in the 
sector above |vaci^_^) (see the discussion around (|3.9p ). That is, if a; = 0, it is valid for excited 

states in the NS vacuum, hence with nj,nj G Z + ^. For = ^, it is valid for excited states in 
the R vacuum, hence with nj, fii G TL. 

When O is a twist field, its associated branch cut changes the sector of the bra or the ket, 
hence formula (|7.18p can then be applied only if one of the bra or the ket is the vacuum, and 
if the branch cut associated to the twist field is chosen so that this vacuum is in the opposite 
sector (in order to keep the excited states in the same sector). If u; = 0, the vacuum will then 
be in the R sector, and if a; = |, it will be in the NS sector. For a branch cut to the right, it is 
the bra that must be chosen as this vacuum, whereas for a branch cut to the left, it is the ket. 

It is easy to check, using (I7.18p . that the formulas above for finite-temperature form factors 
reproduce the known form factors on the circle [191 [20l [2T] . 

7.6 Two-point functions, Predholm determinants 
and scaling limit of the quantum Ising model 

The finite-temperature form factor expansion (|5.5p now gives explicit expansions for finite- 
temperature two-point functions of twist fields at x > 0: 

((cj+(a;,t)cj_(0,0)))^ = e(^[H-£[^H)^ 

k 

Ej (ima: sinh 9j— jmt cosh 6j ) 

d^i • • • dek&=^ 



X 

k 



j=l 

'ilrF^H0r,L))' n tanh(t^)'"^ (7.19) 



and 



((M+(x,t)^_(0,0)))^ = -e(^H-£[|+-]> 



oo „ 



X " 



t j (imx siah dj—imt cosh 9 j) 
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xi'UrFi;-{efL)f n tanh(^)'"\ (7.20) 

j=l l<i<j<k ^ ^ 

where a; = or a; = |, and we recall that £[0] = £^r, and ^[1/2] = ^ns are given in (I7.7p . In order 
to fully clarify the meaning of these finite-temperature correlation functions, we recall also that 
at imaginary time t = ix and at positive x = r, they correspond to the following correlation 
functions in the quantization on the circle: 

(((T+(r,ix)cr_(0,0)))^ = /3(vac^|cr(x,r)cr(0,0)|vac^)/3, 
((//+(T,ix)^_(0,0)))^ = /3(vac^^|^(x,r)^(0,0)|vac^)/3, 

where the vacuum is in the R sector if u; = 0, and in the NS sector if w = 1/2. 

Following [2] , where techniques from [461 [17] were borrowed, Fredholm determinant represen- 
tations can now easily be obtained for two-point functions from the formulas 



Uj — u 



det 

i,j I Ui + Uj 



TT ( — — ^ ) k even, 

+ (7-21) 

k odd 



and 



det,J^— =^-^ TT . (7.22) 

Formula (j7.2ip gives 

((cT+(x,t)a_(0,0)))^ = det(l + K), 

where K is an integral operator with an additional index structure, defined by its action 
(KfUe) = E JZodO'K,^e'{0,e')fAO') and its kernel 

e' = ± 

3/ n\ pelimx s'mh 9 —imt cosh 6) 

K,^,> {e, 9') = i{^F^^+ {9; P)Ytanh' ^ 



^ _j_ g27riajg— em/3 cosh 



Finally, in order to obtain two-point functions of disorder fields, we must consider the linear 
combinations a ziz fi. Formula (j7.22p gives 

(((a+(x,t) +7?/i+(x,t))(a_(0,0) +r?/i„(0,0))))^ = det(l + j(^)) 
with r] = ± and by definition {J^^^f)^{u) = J2 fo° du'J^^J,{u,u')f^'{u') where the kernel is 

e'=± 

given by 



1 |(im/r(u— u ^)—imt{u+u ^)) 



1 + 



The interest in Fredholm determinant representations is, in part, that they can be used to 
efficiently obtain asymptotics of correlation functions. 

Finally, we mention that these two-point functions in the Majorana theory can be used to 
evaluate the off-critical scaling limit of two-point functions in the quantum Ising chain (see, 
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for instance, the book [H]). The quantum Ising chain is a quantum mechanical model with 
Hamiltonian 

-falsing = - ^{JSjSj+i + hs^'j) 

i 

with J > 0. The spin variables s| and are in the spin-1/2 representation of SU{2), and are 
two of the usual Pauli matrices on the j^^ two-dimensional space, the third one being s^: 

10^' J' \o -I 

It is the "Hamiltonian limit" of the two-dimensional Ising classical statistical model. There is 
a value h = he of the transverse magnetic field at which this model is critical. The conformal 
field theory that describes it is the free massless Majorana theory. For h < he, the system is 
ordered, and at zero temperature the average of Sj is non-zero. On the other hand, for h > he, 
the system is disordered. As h is made to approach he, the correlation length ^ associated to 
the two-point function {sjSq) diverges. The scaling limit is obtained by looking at the situation 
where h — > he, while the inverse temperature is made to diverge as J/3 oc and the distances 
between points in correlation functions are made to diverge as \j\ oc The quantum field theory 
model that describes the appropriately normalised correlation functions obtained in this limit is 
the free massive Majorana theory, the product of mass times position being equal to mx = 
If h is sent to he from below (ordered regime), then we have the correspondence 



where 



and Z is a non-zero, non-universal number. On the other hand, if h is sent to he from above 
(disordered regime), then 



Z 



-1 



J Tr (e-^^M"gs|(t)sg(0) 



TV (e-/^^r^in.) 4((;,+ (x,t)/x„(0,0)))|. 



It is important to realise that the spin variables does not converge, in the scaling limit, to 
the twist fields a, fx; indeed, only its products converge to products of twist fields. This is clear, 
since the finite-temperature average of single twist fields are non-zero (but have non-trivial space 
dependence, as explained in Subsection 13. Sp . but finite-temperature averages of spin variables 
are zero (since at finite temperature, there can be no symmetry breaking). One should recall 
that the passage from the quantum Ising model to the Majorana theory involves writing the 
spin variables as exponentials of sums of (bilinear of) fermionic variables lying on a segment of 
the chain, and the two end-points of the segment correspond to two spin variables. 



8 Perspectives 

We have developed partly the concept of finite-temperature form factor in the general context of 
factorised scattering theory, and we completed the program in the case of the Majorana theory. 
The most important next step is, of course, to complete this program in models with non-trivial 
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scattering. We believe that ideas concerning the relation between finite-temperature form factors 
and matrix elements in the quantization on the circle will lead to restrictive conditions that will 
greatly help fully fix finite-temperature form factors in interacting models. Also, the operator 
implementing the generalisation of CFT's "mapping to the cylinder" may be useful, and this 
method is not far from the explicit construction of "boundary-creating operator" in integrable 
boundary QFT. The generalisation to interacting models is a very important step, and will 
open the way to results about large-distance and large-time behaviours of correlation functions 
in interacting, integrable models. 

Another interesting avenue is to generalise the program to the free Dirac theory; this should 
not pose any difficulties, and will clarify the structure of finite-temperature form factors of more 
general twist fields (two-point functions at finite-temperature were already studied in |49| ) . 
Then, it would be interesting to understand the structure for descendants of twist fields in such 
free fermionic models, perhaps using the operator Q defined in (j6.7p that provides a "mapping 
to the cylinder". 

Finally, one would like to obtain the full large-time expansion of correlation functions in the 
quantum Ising model. Besides directly using the finite-temperature form factor expansion, it is 
possible that the Fredholm determinant representations obtained here can be used fruitfully for 
this purpose (work is in progress |31j). 

A OPE's in the Majorana theory 

The order and disorder fields have operator representations a± and fi± on 7i. These operators 
are completely characterised by the leading terms in their OPEs with the fermion fields: 

i;ix,t)a+iO,t) ~ , , ^ u+iO,t), a+{0,t)4:{x,t) ~ , \ »+(0,t), 

2v— vrx — zO+ 2y/—7rx + i0+ 

^l;{x,t)a.{0,t) ~ —-J===^,_{0,t), cT_(0,t)V'(x,t) ~ —-J===f,_{0,t), 

ZVTTX + ^0"•" IVTTX — lU~^ 

iPix, t)fi+{0, t) ~ -—^L==a+{0, t), /i+(0, t)i;{x, t) ~ -—J===a+{0, t), 



V^(x,t)^_(0,t) - , \ a^iO,t), fi^{0,t)i;{x,t) ~ -_^L==a„(0,t) 



and 



i;{x, t)a+{0, t) ~ -—^==^+(0, t), a+{0, t)i^[x, t) ~ -—^==^+(0, t), 
2v— vrx + zO+ 2v— TTX — zO+ 

V;(x,t)a_(0,t) ~ - ' /»-(0,t), a_(0,t)V5(x,t) ~ --_^==^_(0, t), 
2v vrx — zO+ 2v vrx + zO+ 

i^{x, t)fi+{0, t) ~ —=L==a+{0, t), fi+{0, t)ij{x, t) ~ --^L==a+{0, t), 
2v— 7rx-|-i0+ 2\/ —TTX — tO+ 

^P{x,t)fi.iO,t) ~ ^ tT-(0,t), ^L-{0,t)ij{x,t) - -—L=a_{0,t), 
2v7rx — iO~*~ 2\'TTx + 

where everywhere, the square root is on its principal branch. 
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